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Abstract. 

The  problem  of  quantization  of  the  electromagnetic  field  in¬ 
side^  resonator  filled  with  a  dielectric  medium  with  time- 
dependent  characteristics  in  the  presence  of  4;he ^external  alter¬ 
nating  current  is  studied.  The  exact  propagator,  coherent  and 
Fock's  states  are  obtained  for  a  separate  quantized  mode  of  the 
field  represented  by  a  quantum  oscillator  with  time-dependent 
frequency  linearly  coupled  via  momentum  with/, external  current. 
Delta-pulse  excitation  of  the  medium  is  considered  as^^exampie. 


1.  Introduction. 


Time -dependent  systems  are  observed  in  various  physics  expe¬ 
riments.  Two  general  types  of  such  systems  are:  that  which  is 
formed  through  its  own  environmental  conditions,  and  that  which 
is  formed  when  external  forces  are  added.  In  regard  to  the  second 
type,  various  experiments  are  being  carried  out  to  see  how  an  ap¬ 
plied  time -dependent  electric,  magnetic  or  other  field  can  alter 
the  physical  properties  of  materials  such  as  semiconductors  and 
superconductors.  Experiments  show  that  a  system  becomes  time- 
dependent  when  a  time-dependent  electric  or  magnetic  field  Csuch 
as  a. c. ]  is  appi led. 

The  aim  of  the  present  paper  is  to  investigate  the  behaviour 
of  a  quantized  mode  of  the  electromagnetic  field  inside  a  resona¬ 
tor  filled  with  a  dielectric  medium  acted  on  by  some  external 
"pumping"  electromagnetic  field  and  external  alternating  current. 
Due  to  this  "pumping"  the  properties  of  the  internal  medium  beco¬ 
me  time-dependent.  ^ 

First  we  study  the  problem  of  quantization  of/^electromagne- 
tic  field  in  the  case  of  a  time-dependent  medium  and  external 
current.  We  show  that  the  dynamics  of  every  mode  is  the  same  as 
the  dynamics  of  a  time-dependent  forced  harmonic  oscillator.  There 
4Wlfore,the  rest  of  the  paper  is  devoted  to  the  problem  of  finding 
exact  solutions  for  the  time -dependent  forced  harmonic  oscilla¬ 
tor.  This  last  problem  was  considered  by  many  authors;  see.  e.g. . 
refs./l-13/  and  especially  /14-19/  where  quantum  systems  with  the 
most  general  quadratic  multidimensional  Hamiltonians  were  studied 
in  detail.  The  difference  of  our  study  from  the  previous  ones 
consists  in  the  unusual  coupling  with  the  external  force;  not 
through  coordinate  but  via  momentum.  Besides,  we  represent  the 
known  formulas  for  the  propagators  and  eigenfunctions  in  another 
parametr ization  than  in  the  previous  studies  and  calculate  some 
products  of  matrix  elements  which  were  not  given  earlier. 

The  special  case  of  time-dependences  of  the  frequency  and 
external  force  in  the  form  of  delta-pulses  is  considered  to  il¬ 
lustrate  general  formulas. 


3 


2.  Field  quantization  in  a  time-dependent  medium  in  the 
presence  of  an  external  current. 


■i/,  I 

The  well-known  usual  scheme  oi  quantization  of/^eiectyqr/raagne- 
tic  field  in  trhe^uniform  time- independent  medium  Csee.  e.g./20/D 
results  in  the  replacement  of  the  field  variables  by  a  countable 
set  of  generalized  canonically  conjugate^  coordinate  and  momentum 
operators  obeying  harmonic  oscillator  equation  of  motion.  Here  we 
consider  a  more  generai  case  when  dieiectric  and  magnetic  permea¬ 
bilities  depend  on  time  Cbut  not  depend  on  space  variables^,  and 
some  time-dependent  externa i  current  is  present, .  The  same  prob¬ 
lem.  but  without  externai  current^  was  considerea  in  /21-'.  and  the 
most  general  case  of  nonuniform^  and  time -dependent  medium  C  al¬ 
so  without  current)  was  studied  in  /22/. 

For  the  sake  of  simplicity  we  confine  ourselves  to  the  simp¬ 
lest  case  of  the  so-called  "unidimensionai  electrodynamics",  when 
linearly-polarized  mutually  perpendicular  electric  and  magnetic 
fields  depend  only  on  the  single  space  variable  ^  and  time  t. 

Then  Maxwell’s  equations  assume  the  following  form: 
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Introducing  the  vector  potential  as  usual . 
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we  get  the  second  order  equation 
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so  that  the  Ham.itonian  density  is  given  bv 
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We  assume  the  field  to  occupy  an  ideal  resonator  with  walls  at 
the  points  ^=0  and  ^=L.  Then  the  following  mode  decomposition  is 
natural : 
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After  this  the  field  Hamiltonian  turns  into  the  Hamiltonian  of  a 
countable  set  of  noninteracting  harmonic  time-dependent  forced 
osci 1 lators; 
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To  auantize^ami  Roman  C2.  lOD,  we  have  to  treat *var lables  x 
and  p^  as  operators  satisfying  canonical  commutation  relations. 
At  least  two  variants  are  possible.  First,  we  may  consider  p^  as 
the  aeneralized  coordinate,  and  x  as  a  Generalized  momentum.  Se- 

I) 

condly.  we  mav  make  the  opposite  choice:  as  a  generalized  mo¬ 

mentum.  and  X  as  a  Generalized  coordinate.  In  both  cases  we  have 

n 

the  harmonic  oscillator  Hamiltonians  v,'ith.  generally  speaking, 
variable  masses  and  frequencies.  Now  let  us  take  into  account 
that  really  the  magnetic  permeability  differs  from  unity  to  a  very 
small  value  in  all  cases  when"^  inear  equations  C2.  D  hold. 
Thus  from  the  physical  point  of  view  it  is  sufficient  to  consider 
the  function  ^iCt)=const=l.  Then  we  mav  study  either  the  oscilla- 


tor  with  a  constant  frequency  but  with  a  time-dependent  mass,  or 
the  oscillator  with  a  constant  mass  and  with  a  time-dependent 
frequency.  The  second  situation  seems  much  more  familiar  and  con¬ 
venient.  Therefore  our  choice  is  to  treat  as  the  generalized 
coordinate  operator,  and  p^  as  the  generalized  momentum  operator 
satisfying  the  commutation  relations 


C2. 12) 


The  physical  significance  of  such  a  choice  is  that  the  electric 
displacement  vector  plays  the  role  of  the  generalized  coordinate 
while  the  generalized  momentum  is  related  to  the  vector  potential 
or  the  magnetic  induction  vector  (see  also  refs.  .''23-26/  in  this 
connection). 

Since  ail  .modes  in AHami Itonian  C2. 10)  are  uncoupled,  we 
shall  omit  hereafter  the  number  of  the  mode  and  coTnsider  the 
following  one-mode  Hamiltonian;  ^ 


[  p^  +  (/ct)x^  j  +  rct)p  . 


C2. 13) 


_  CcTm)^ 

t/ct)  =— - - 

L^cCt) 


C2. 14) 


This  Hamiltonian  possesses  the‘^cur icus  peculiarity;  it  contains 
coupling  with  an  external  time -dependent  current  C  represented  by 
function  fCt)  )  not  through  the  coordinate,  but  through  the  mo¬ 
mentum.  Therefore, obtaining  explicit  solutions  to  this  not  very 
usual  oroblem  seems'^rather  interestinq  task. 

A 


3.  Propagator. 


We  begin  with  calculating  the  exact  propagator  for  the 
Schrodinger  equation,  i.e.  the  integral  kernel  relating  the  valu¬ 
es  of  the  wave  function  at  the  instants  of  time  t’  and  t 

^^rCx.t)  =  r  dx'KCx.t;x't’)^Kx’ .t’) .  t>t'.  C3.1) 

Since ^Hami Roman  C2. 13)  is  a  nonuniform  quadratic  form  of  the 
coordinate  and  momentum  operator,  the  propagator  also  is  an  expo- 


nentiai  of  some  nonuniform  quadratic  form  --12-19/ 

KCx.t;x’.t’)  =  exp^aCt.t')x^  +  bCt.t'Dxx'  +  cCl.l’Dx’^  + 

+  gCt.t’Dx  +  hCt.t'Dx'  -  dCt.t')  j.  C3.2] 

Substituting  this  expression  into  the  Schrodinger  equation,  one 
can  obtain  a  system  of  coupled  nonlinear  ordinary  differential 
equations  for  time-dependent  coefficients  of  the  quadratic  form, 
e.g. 


=  2Lha^  -  w^CtD  /2Lh  . 


C3.3] 


—  =  2Lfiab. 
dt 

do  ih  _ 

—  =  —  b^. 
dt  2 

do  i 

—  =  -  2ihc^  ^  —  oj^Ct’D. 

dt'  2h 
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2Lhcb, 
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r  o 

V  O  .  - 


and  similar  equations  tor  the  coefficients  of  lin^.r  form  in 

C3.21.  Eqs.  C3.33-C3.S1  can  be  solved  with  the  aid  of^anzatz 

L  dq  i  da 

act)  = - ,  cCt'3  =  - - -  .  C3.93 

^  2hqCt)  dt  2hqCt')  dt’ 

where ^new  function  qCt)  obeys  the  classical  equation  for  the  os¬ 
cillator  with  a  T  ime-deoendent  ■''requencv 


d"q  . 

—  +  ofCtlaCt'  =  0. 

Z 

U 


C3.  10) 


However,  certain  care  is  requi^d.  since  eq.  C3. 10)  possesses  two 
independent  solutions,  so  t  ha  t^i  unci  ions  aCt)  and  qCt’.'  in  C3.9j 
are  different  in  the  aeneral  case. 

There  exists  a  more  direct  and  convenient  method  of  calcuia- 
tinq  the  propagator  oroDosed  in  /16/  and  developed  in  /IT.IS.'O  It 


IS  based  on  the  concept  of  time -dependent  integrals  of  motion 
Cwhose  role  for  obtaining  exact  solutions  of  quantum  mechanical 
problems  was  stressed  in  refs. --'1 .2/3 .  An  integral  of  motion  is  an 

operator  satisfvina  the  equation 

*  ^ 

31  p  n 

ih  —  =  I  H.I  .  C3.113 

3t  >-  J 

Suppose  two^operators  XCt.t']  and  PCt.t’3  satisfying  the  conditions 
i{Ct.t’3=x,  PCt.t’}=p  are  known.  Then  the  kernel  KCx.tix’f)  can 
be  found  from  the  equations  -'IR-IS.' 


XCt.t’)KCx.t;x’t’)  =  x’KCx.tix't'). 


PCt.t’}KCx.t.x’f:  =  ih—  iiCx.tix’t 


3 

1— 
3x’ 


:3.  123 
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Pryr  quadratic  Hamiltonians  ooerators  X  and  P  are  evidently  linear 
combinations  of  operators  x  and  p: 


PCt.t’3  =  ;\^Ct.t'3p  +  AgCt.t'Dx  +  5^Ct3, 

XCt.t’3  =  A^Ct.t’Dp  +  A^Ct.t’Dx  +  S^Ct). 

Substituting  these  exoression  into  C3. 113  we  get  a  set  of 
linear  ordinary  differential  equations 
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—  =  -  A„.  — ^  =  u^CtlA,, 

dt  2  dt  ^ 

dA  dA 
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dt  ^  dt  ^ 

with  the  initial  conditions 
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L  163 


C3.  173 


A^Ct.tJ  =  A^Ct.tj  =  1.  AgCt.t;  =  A^Ct.tl  =  0. 

6^Ct.t3  = 

The  solutions  of  eqs.  C3. 163-C3. 133  are  as  follows; 

dq,  \  UM 

A.Ct.t’:^  =  q.ct.f'.  A  rt.f'  -  -  — ^  5At.f  =  ffCTj  — 

^  ^  ,  dT 


dq, 


dx. 


C3. 193 


^  5,Ct.t’)  -  -  ffCT]— ^  dz. 

^  ^  dl  ^  {.dr 

C3.20} 

where  both  functions  and  satisfy  the  same  equation  C3.10D 
but  with  different  initial  conditions 

q.Ct.tJ  =  1.  — V.t.t)  =  0.  q_Ct.t)  =  0.  — ^Ct.tD  =  1.  C3.21.') 

^  dt  ^  dt 

After  this, equations  C3.12D  and  C3. 133  immediately^  determine  all 
coefficients  of  the  quadratic  form  C3.23  except?feg(|jj^j^the  free  term 
fCt.t’;  which  can  ne  found  from  rhe  Schrddinaer  equation 
dK  .. 

ih—  =  HK.  C3.22D 

dt 

The  final  result  is  as  follows; 

KCx.tix’.t’)  =  ^27rLhq2Ct.t'3j  exp - -  [  2xx']- 

C3.233 

■  K  ^  s  J  fCT)qc.5F,CT) 


wnere 


FjCt.t’3  =  JfCTjqjCrldT,  j=1.2. 


Now  let  us  discuss  the  structure  of^functions  q^  and  q^  -- 
solutions  of  eq.C3. 103.  Suppose  we  look  for  a  complex  solution  of 
this  equation  in  the  form  ^  qCt3=T7Ct3exp[  iyCt3  ] ,  77  and  y  being 
real  functions.  Then  we  aet  two  eouations 


•  •  p  •  p 

7)  +  Q'^CtJT)  =  TTy*^. 

w  +  27^=0 

From  the  last  equation  we  get 
rfr  =  const  =  n. 

Then  fifct ion  TjCtJ  obevs  the  nonlinear  equation 


C3.253 

C3.263 


C3.273 
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T)  +  w^ct)7?  =  nP/7)^.  C3.28D 

Suppose  T?Ct3  is  an  arbitrary  solution  of  eq.C3.283.  Then „ function 
«35es^  q^Ct. t '  3  and  q^Ct.t'Dcan  be  expressed  in  terms  cf  7)Ct3  and 

ndT 

yCt.f)  =  I  -  C3.29D 

as  follows: 
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Using  these  expressions,  we  can  rewrite  C3.233  in  a  more  symmet¬ 
ric  form: 
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4.  Coherent  states. 

If  we  choose ^comDle:-:  solution  of  eq.  C3. 10)  satisfying  an  ad¬ 
ditional  condition 

qq»^  -  q^q  =  2i,  C4.  1) 

which  IS  equivalent  to  the  choice  £2=1  in  relation  C3.27).  then  we 
can  construct  a  nonhermitian  operator  --linear  integral  of  motion 
n4.i8/ 


ACt3  =  iCZhD  qCt^p  -  qCtDx  +  JdTfLTDqCT) j  ■  =  C4.2D 

t 

=  |^7Kt::p  -  Ch  +  lyTjJx  -  J  drfCT^LT?  + 

t  ’ 

-  ]  U.S 


:4.3i 


Cl'  IS  arbitrary)  satisfying^ bosonic  annihilation-creation  opera¬ 
tors  commutation  relation 


[  ACt) , 


A  Ct)  =  1 


C4.4) 


Thus  we  can  ccnstruci  an  overcomplet-='  set  cf  conerent.  states 
iir  Cx.t)  satisfying  both  the  Schrddinger  equation  and  the  equation 


ACD^A^Cx.t)  =  o^lr^Cx.t) . 


C4.5) 


a  being  an  arbitrary  complex  function.  The  explicit  expression 
for  these  states  can  be  written  as  a  special  case  of  general  for¬ 
mulas  for  coherent  states  of  multidimensional  quadratic  systems 
found  in  /15.13/ 

r  r  (  7  {  ih  \  ^  I 


y  -1I/4  f  y  f  ]  y  1 

—  exD-  -  -  1  -  —  X  -  -  lyCt. 

nh  J  ■  1  2h  i  w  }  2 


1  r  y  r''"  r 

1  |a|2  .  IJ  VSK-a  - 

2  I  f>h  L 


'  d  <•  1  1 

ifdxfCT)  —  I —  "■  + 

f  ^^1  ''IcTT  J  1 


2  1’^  ,»(*!  r  d  fslnjrCt.Ti 

hj  f 


^  r 

-  I  dxfC 

J 


d  fsinyCt.T)!  I  d  f  e'^yct.s) 


yCx) 


fdsfCs)  — 

J.  ds 


The  expectation  values  of  x  and  p  in  these  states  are 

f  2h  d  fsinyCt.x) 

<a|x|a>=  — Re  -  f  dxfCx) — - : - : - : 

[  y  J  ^  ^  dx  [  [yCDyCx)] 


.C4.7) 


C4.83 


d 

<a|p|a>  =  —  <a|x |a>. 

dt 

The  variances 

(7p  =  <a|p^|a>  -  ^ca|p|a>j  .  =  va|x^|a> 


aepena  neither  on  a  nor  on  fCt).  but  are  compieteiy  determined  by 
the  function  qCt]  /18/: 


X 


p 


IqCt 


"I  = 


h 


^  •  2  ^  -r 

-  IqCt}  1^  =  -  r [  i 
Their  product 

] 


C4.9) 

C4.  lOD 


C4.1D 


vanes  in  time;  moreover,  it  is  greater  than  the  minimal  possible 
value  h‘'/4  provided  ifK).  This  is  explained  by  the  fact  that^state 
(4.6)  IS  really  correlated  coherent  states  /18.27/  with  nonzero 


correlated  coefficient 


<xp  +  dx>/2  -vx><p> 


O'  a  1 
X  p 


1/2 


.  •  2  ^ 
If)  +  r 


^  1/2 


''-4. 123 


Correlated  coherent  states  minimize  the  Schrodinger -Robertson  gen- 
^eraiized  uncertainty  relation  /13.27/ 


1  ^  C4.133 

Since  the  substitution  of  C4.11J  and  C4.12j  into  C4.133  trans¬ 
forms  this  inequality  into  the  identity. 


5.  Fock's  states. 


It  is  well  known  /23/  that  coherent  states  are  generating 
functions  for  Fock's  states: 


12 


1“^  “  n>. 

n=0  C  n !  ) 

C5. 1) 

A  A 

A  A|n>  =  n  |n>. 

C\J 

tri 

On  the  other  hand,  on  exponential  function  of  a  quadratic 
form  IS  the  generating  function  for  Hermite's  polynomials  /29/. 
Thus  one  can  obtain  the  following  formula: 

•  x|n>  =  ^^Cx.tD  =  [s'Vilj 


o  H 

n 


f  dTfCr) — 

.  dx 


sinyCt.T^ 


yCtDyCxD 


-V 1/2 


C5.3D 


^^^Cx.t)  being  given  by  eq.C4.6)  with  a=0.  Functions  C5.3)  are 
also  eigenstates  of  the  operator  (considered  first  by 


Lewis  /I. 2/) 

ECtJ  =  hnj^A^A  +  1/2]  =  ^  -  [ 

^  /V  /\  A. 

-  7J7?(XP  +  pxj  + 

2Re^FCtje~^^^^^  |^77^p  -  lyxj  +  7)xj 

where 

t 

]  *  |FC13|^' 

|-  .  C5.4) 

FCt)  =  J  dxfCxl^T?  +  lyTjje^^rCT)^ 

C5.5) 

and^functions  iKtli.  yCt)  satisfy  eqs.  C3. 25D-C3. 29D  with  the 
given  positive  constant  f2. 


ECt)<^^Cx.tj  =  hn[n  ^  •  n=0.1.2....  C5.6^ 

Tkx  Operator  ECt)  is  an  integral  ofy^motion  coinciding  with  the 
usual  energy  operator  in  the  case  of  time- independent  functions  £2 
and  f  C  provided  the  choice  £2=w  is  made).  Therefore.it  can  be 
cal  led /^"general  1  zed  energy  operator". 

If  we  define  an  "uncertainty  product"  as 


d.3 


AxAd  = 


p  “  {[[ 


V  ^  I  n  >  -  < 


<m  X  n> 


m 


|x|n>^  j  ^  <m|x^|n>  -  <m|xln>^  jj 


<m|p^|n>  -  <m|p|n>^  j  <m|p^|n>  -  <m|p|n>^  jj  j-  , 
then  the  following  formulas  can  be  obtained: 


Hp]->.n=  K 


■rf  ^irzr  . 

1  -  J  [n  +  I/2J 


•  2  2 
y  T 


hH-2.n=  '■[l  '  "  2Kn  .  n] 

9  U 

All  other  products  are  equal  to  zero. 


1/2 


1/2 


5.  Mode  excitation  by  delta-impulse. 


1/2 

o 


C5.7D 


C5.8D 

C5.9) 

C5.  ICO 


As  an  example, we  consider  the  simplest  time  dependences  of 
the  effective  frequency  and  "force"  in  the  form  of  delta- impulses: 


u^Ctj  =  ^  -  mi:.  W>0. 


fCt)  =  ] 


f 


|t|  1  a/2 


L  0,  |t|  ;■  a/2 


a^.  f-x®,  af-»F=const. 


C6. 1.' 
C6.2) 


The  dependence  C6. 13  can  be  realized  approximately  in  the 
case  of  ^  very  fast  ionization  of  the  medium  when  the  dielectric 
permeability  e  =  1  -  47me^/mu^  turns  into  zero  due  to  an  instant 
growth  of  the  electron  concentration  l  here  10  is  the  resonance 
frequency  of  the  mode)  and  then  quickly  restores  -ihi?-)finitial  value 
due  to  recombination.  Some  physical  effects  in  such  a  "plasma 
window"  were  discussed<iS5(#Si^  in  ref .  ■'^O./. 

Further  calculations  are  based  in  part  on  results  of  ref. 
/31./^ 

FThe  Hami It  Ionian  of  this  system  is  of  the  form 


// 


H  =  _ —  +  -  +  W5Ct]  +  rCt}p.  C6.3D 

2m  2  ^  J 

— -i^oi  lowing  the  usual  procedure  /iS-''  and  taking  into  conside¬ 
ration  the  results  obtained  in  /?A/  onp  can  construct  the  linear 

TJd- 

integral  of^ mot  ion  for  the  Hamiltonian  C6.3D. 


/V 


sinOt 


iOe 


L£2t 


-  Wcosnt 


U  6CtD. 


t>0. 

where 


C6.4D 


p^  -  x^=  [h/mn]^''?  5Ct>[-^YinF  -  FW/2] .  F  =  limf^a. 

f  -XD 
O 

-ui. 

One  can  checkin  that^  integral  of  ^motion  C6.4)  and  its  hermi- 
tian  conjugate  operator  satisfy  the  commutation  relations  of  bo¬ 
son  creation  and  annihilation  operators, 


[A. A*]  =  1. 

a 

The  ground  state  of^^ quantum  oscillator  transforms  to  a  cor¬ 
related  coherent  state  C4.63  after  5-kickin^of^^reauency  C:t  is 
necessary  to  out  aCt)  =  e*"^  -  CW/£23sinJ2t  in, formula  C4.23  3.  with 
the  variances  C4.93.  C4.10)  equal  to 


^  r 

2mnt  ^ 

%  =  — ( 1 

P  2mn^ 


-„s  1  n^J2t 
-  pos^nt 


w 

Q 

W 

a 


sin2ntj . 
s i n2nt  j . 


The  correlation  coefficient  C4.123  is  nonzero  after  5-kicking 
and  IS  equal  to 

r  =  i  -  I  1  +  [w^4n^]sin22nt  +  [w^jcos^£2t  +  [w?2nP j  s  i  n4nt]-  \ 

The  squeezing  i^fficient  is  not  equal  to 

unity  and  is  given  bvv'ormuia 


k  = 


T.r2 

**  ? 

1  "  -  sin  J2L 

^  T" 

1  +  -_cos  nt 


1/2 


-  3in2£2L 


w 

-  sin2£2t 
Q 
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